3->"»fc 




RESEARCH DEPARTMENT 



Calculation of the field strength required 

for a television service, in the presence 

of co -channel interfering signals 

Part 3: The computer programme 
RESEARCH REPORT No. RA-12/3 

UDC 621. 391. 812. 8. 1968/46 

681. 3.06 



THE BRITISH BROADCASTING CORPORATION 
ENGINEERING DIVISION 



RESEARCH DEPARTMENT 

CALCULATION OF THE FIELD STRENGTH REQUIRED FOR A TELEVISION SERVICE, 
IN THE PRESENCE OF CO-CHANNEL INTERFERING SIGNALS 

Part 3: The Computer Programme 



Research Report No. RA-12/3 
UDC 621.391.812.8 1968/46 
681.3.06 



This Report is the property of the British 
Broadcasting Corporation and may not be 
reproduced in any form without the written 
permission of the Corporation. 

It uses SI units in accordance with B.S. 
document PD 5686. 



R.W. Lee, M.Sc. Head of Research and Development 



Research Report No. RA-12/3 
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CO-CHANNEL INTERFERING SIGNALS 

Part 3: The Computer Programme 



PREFACE 

The report on 'Calculation of the Field Strength Required for a 
Television Service, in the Presence of Co-channel Interfering Signals' 
is in three parts, as follows: 

Part 1: The assessment of a single interfering source 

Part 2: Effect of multiple interfering sources 

Part 3: The computer programme 

Parts 1 and 2 deal with the theoretical background of the sub- 
ject, and have been written as entities which can be read separately; 
for a full understanding of the work, however, the reader is advised to 
study Parts 1 and 2 together. Part 3 describes the form and flow diagram 
of the original computer programme to carry out thecalculation processes 
described in Parts 1 and 2; recent improvements in detail and adapta- 
tions to a larger and faster computer have been made, but are not con- 
sidered here. 
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CALCULATION OF THE FIELD STRENGTH REQUIRED FOR A TELEVISION SERVICE, 
IN THE PRESENCE OF CO-CHANNEL INTERFERING SIGNALS 

Part 3: The Computer Programme 



1. GENERAL 



SUMMARY 

Parts 1 and 2, which have been issued separately, give the theory 
underlying the calculation processes which have been actually used in the 
BBC Research Department. Part 3, which follows, describes the original 
computer programme used for carrying them out. 



(viii) The four components resulting from (iv), (v), 
(vi) and (vii) are summed to produce the 
protected field strength. 

(ix) The foregoing procedure is repeated for all 
unwanted sources. 

(x) With the medians and standard deviations 
derived in (ix) a combined time-location 
distribution is produced for each unwanted 
source. 

(xi) These are combined by multiplying prob- 
abilities at chosen field-strength levels. 

(xii) The resulting distribution is divided into 
segments each approximating to a Normal 
djstribution; the median and standard deviation 
in time-location for each segment are deter- 
mined. 

(xiii) For each of these segments a probability at 
the mi d-fi eld-strength level is derived by 
multiplying the probabilities due to location 
variation only. 

(xiv) From (xii) and (xiii) the standard deviation in 
location variation only is calculated and the 
standard deviation in time deduced. 

(xv) These are used to derive the 95% time 50% 
location, 95% time 70% location protected 
fields and the time percentage for which 70% 
location will be protected for field strengths 
of 70 dB and 80 dB(pV/m). 



Before describing the computer programme in 
detail, it is useful to state the stages of the cal- 
culation necessary for the determination of co- 
channel interference levels, and hence the protected 
field strength in the presence of multiple interfering 
sources. 

For each receiving test point we first calculate 
the interference due to individual co-channel inter- 
fering signals, as described in Part 1 1 of this report. 
We then calculate the protected field strength in the 
presence of a single interfering signal, and finally 
the protected field strength in the presence of all 
the individual interfering signals taken together, as 
explained in Part 2 2 of this report. Throughout the 
report, field strengths and standard deviations are 
regarded as expressed in decibels. 

The procedure is as follows: 

(i) The receiving location to be examined is 
selected. 

(ii) Co-channel sources are listed. Then, for 
each of these: 

(iii) The distance to the test point and the bear- 
ings to and from the test point are found. 

(iv) The carrier -off set relationship of the wanted 
and unwanted signals is derived and the 
resulting protection ratio found. 

(v) The bearing from the test point to the wanted 
transmitter is compared with (iii) to produce 
a receiving aerial directivity allowance. This 
is combined with a polarization allowance (if 
any). 

(vi) The e.r.p. of the unwanted source (in dB rela- 
tive to 1 kW) is obtained. 

(vii) The median field strength of the unwanted 
signal for an e.r.p. of 1kW is calculated, 
together with the time variation. 



2. THE COMPUTER 

The computer originally used for this programme 
was an Elliott 803B digital computer with a 8192- 
word core store, an automatic floating point unit and 
two film handlers. 

Input and output are both effected by means of 
punched paper tape. 



The films used are 35mm magnetic films with a 
storage capacity of 4096 64-word blocks on each 
film. Each computer word has a fixed length of 39 
binary digits (bits). 

Recently the number of transmitters to be 
considered has been increasing rapidly, and the 
Elliott 803B computer is no longer adequate. The 
programme has therefore been adapted for an ICT 
1905 — a much larger and faster machine with a 
32,000-word store and eight magnetic-tape decks, 
having paper-tape input and output to a line printer. 
Considerable improvements in the detail of the 
calculation have been achieved, but these are not 
discussed here, since the principles used and objec- 
tives sought remain unchanged. 



One film (A) contains five separate types of 
information: 

(i) Data for the coast-line grid 
(ii) A decode table for (i) 
(iii) Transmitter/ channel correspondence list 
(iv) Transmitter details 
(v) Transmitter and country names 

A second film (B) contains the test location 
details. 

A separate programme is provided for the stor- 
age of these items and provision has been made for 
three specific types of change (transmitter location, 
carrier offset and e.r.p.) to transmitter details which 
are already stored. 



3. DATA STORAGE 

It was realized that the amount of data to be 
stored, up to 85 numbers per transmitter and 39 
numbers per test location with eventually some 
hundreds of transmitters and perhaps two thousand 
test locations, was far in excess of that which the 
working store could accommodate. 

It was also felt that input of the required data 
on punched paper tape each time the programme was 
run, was undesirable. Because of this it was 
decided to use films to store the information. 

Since the film block length is only 64 words, 
the information had to be "packed". This involves 
storing two or more numbers in one computer word. 
The transmitter information can then be fitted on to 
one film block as can each set of test location 
information. 



4. FLOW DIAGRAM 

The objective of the programme is to carry out 
the calculations required for obtaining co-channel 
interference levels and hence the protected field 
strength by the methods described in Parts 1 and 
2 . The text which follows includes further details 
of the calculations, where these are required for 
understanding how the programme achieves its 
objective. 

The essence of the required technique is 
explained by means of a flow diagram (shown in 
Table 1). This flow diagram includes references 
to the text of Part 3, where fuller explanation is to 
be found. 

There are two programmes involved in carrying 
out the calculation. The first programme, Table 
1(a), assembles the required data from Films "A* 
and 'B' on a temporary store film. This cuts down 
film search time during the calculation. 



TABLE Ha) 
Data Preparation Programme 



Re). Programme 

1 Coast-line grid, Film A to temporary film store (T.F.S.) 

2 Decode table, Film A to T.F.S. 

3 Country names, Film A to T.F.S. 

4 Read data tape 

5 Channel-transmitter data appropriate to channel specified on data tape, 
Film A to T.F.S. 

6 Transmitter names for channel, Film A to T.F.S. 

7 Repeat, from 5, for all required channels 

8 Sort all transmitters into ascending numerical order 

9 Transmitter data checking for correct channel, Film A to T.F.S. 

10 Repeat, from 9 for all required transmitters 

11 Sort all test locations into ascending order 

12 Replace Film A by Film B 

13 Test location data from Film B 

14 Pack height and store in appropriate transmitter, to T.F.S. 

15 Pack angles and store in appropriate test location, to T.F.S. 

16 Repeat, from 13, for all required test locations 



Reference 
Section 

4.1 
4.2 
4.3 
4.4 

4.5 
4.3 



4.6 



4.7 



TABLE Kb) 



Ref. 

1 
2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 

19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 

30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 

47 



48 
49 
50 
51 
52 
53 



C.C.I. Calculation Programme 

Programme 

Read data tape 

Coast line grid, from temporary film store (T.F.S.) 

Decode table, from T.F.S. 

Country names, from T.F.S. 

Channel -transmitter data appropriate to required channel, from T.F.S. 

Transmitter names for required channel, from T.F.S. 

Transmitter data for service area from T.F.S. 

Check channel; if error — stop 

Unpack test location numbers 

Calculate convergence at transmitter site 

Unpack test location heights 

Terrain angle data for all test locations, from T.F.S. 

Calculate bearings and distance between test location and wanted transmitter 

Calculate path composition between test location and wanted transmitter 

Calculate e.r.p. at appropriate bearing (using linear interpolation) 

Calculate convergence at test location 

Unpack terrain angles at test location 

Calculate terrain angle at test location in direction of transmitter (using linear 

interpolation) 

Calculate terrain angle at transmitter in direction of test location 

Calculate 1% time and 10% time field strengths for the appropriate path 

Calculate 50% time field strength and aj for wanted signal 

Repeat, from 13, for all test locations 

Print offset of wanted transmitter 

Transmitter data for interfering transmitter, from T.F.S. 

Check channel; if error — stop 

Print offset 

Calculate offset allowance 

Unpack terrain angles at transmitter 

Calculate bearings and distance and path composition between test location 

and interfering transmitter 

Calculate convergence at transmitter if in U.K. 

Calculate receiving aerial directivity factor including polarization allowance 

Sum offset and receiving aerial factors 

Calculate e.r.p. 

Sum offset, receiving aerial and e.r.p. factors 

Calculate terrain angles at test location in direction of transmitter 

Calculate terrain angle at transmitter in direction of test location 

Calculate 1% time and 10% time field strengths for the appropriate path 

Calculate 50% time field strength and a T for interfering signal 

Calculate extra protection required for long-term interference 

Derive protected median field strength 

Calculate ay for required signal 

Repeat, from 29, for all test locations 

Repeat, from 24, for all interfering transmitters 

Calculate ct t l for a test location 

Repeat, from 44, for all interfering transmitters 

Calculate 50% location, 95% time protected field strength 

Calculate 70% location, 95% time protected field strength 

Calculate % time protected for 70% locations and 70 dB protected field 

strength. Calculate % time protected for 70% locations and 80dB protected 

field strength 

Calculate 95% time, 50% location protected field for an interfering transmitter 

Repeat, from 48, for all interfering transmitters 

Print out 

Repeat, from 44, for all test locations 

Repeat, from 7, for all wanted transmitters 

Repeat, from 5, for all channels 



Reference 
Section 



4.8 



4.9 

4.10 

4.13 



4.11 
4.12 



.15 



4.17 
4.16 



4.18 



4.18 



Table 2 



The second programme, Table 1(b), describes 
the sequence of calculation showing the transfer of 
data from the temporary store film into the main store 
where required. 

4.1. Coast-line Grid 
The data for the coast-line grid consists of a 
series of distances in kilometres from the North 
Pole along lines of longitude, these lines being 
spaced at 0-1° intervals from 10°E to 10°W. There 
are thus 201 such lines, and the distances stored 



are those corresponding to land/sea intersections 
which occur along each line. To provide a known 
datum, i.e. that the first distance is a crossing from 
land to sea, a dummy coastline at 62°N is included. 
Since the number of intersections along each line 
varies from line to line, a decode table has been 
set up. The total number of intersections to be 
stored is about 1000. These are 'packed' two to a 
word and occupy 9 film blocks, not all of which are 
completely filled. 



TX Site Name: 
N.G.R. 



No. 



IMF Qk and Off. : 



_or Lat. 

L 



"N 



Date: 
Long. 



V.H.F. 



frequencies 
Vision 

Sound 



kc/s 
kc/s 

kc/s 



kc/s 
kc/s 
kc/s 
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Terrain Clearance Angles for TX Ae. 
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310° 




320° 




330° 




340° 




350° 





Rel. F.S. on Bearings East of True North for ('nit Field 

0° 20° 40° 60° 80° 

120° 140° 

260° 



240° 



160°_ 
280° 



180 c 
300° 



200° 
320° 



Test Locations within the Service Area 

No. Name 

1. 



No. 
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340° 



Name 



Fig. 1 • Transmitter data sheet 



4.2. Decode Table 

The decode table gives, for each line, the 
position of the first element of that line in the coast- 
line grid and the number of intersections. These 
two items of information are 'packed' into one 
computer word, so the decode table occupies 201 
locations. For convenience this is stored on film 
and occupies four blocks. 

4.3. Country and Transmitter Names 

Country names are stored in full. Transmitter 
names are abbreviated to a maximum of seven letters 
to facilitate storage in one computer word. 

4.4. Data Tape 

The data tape specifies: 

(a) The number of channels to be considered 

(b) For each channel in (a), the channel number, 
the number of the first United Kingdom trans- 
mitter to be dealt with (treating the transmitters 
on the channel as a set numbered in sequence), 
and the number of United Kingdom transmitters 
to be considered. If required, calculations can 
be carried out on a part of any sequence of 
transmitters occupying the same channel. 

(c) The date 

4.5. Transmitter/Channel Correspondence List 

The correspondence list presently occupies 
one film block per channel. It consists of the num- 
ber of United Kingdom transmitters, the number of 
continental transmitters and a list of these trans- 
mitters in the form of their code numbers. 

This information is stored in a modified form 
in that the transmitter code numbers are replaced by 
the corresponding film block numbers. This is done 
by the storage programme. Channel N occupies film 
block N + 10. 

4.6. Transmitter Data 

Each United Kingdom transmitter has been 
given a pair of numbers. The first is the number of 
the main station and its service area, and this is 
at present in the range 101 — 164. The second is 
for the main station and the relay stations inthe area 
are then numbered in sequence. Provision has been 
made for 20 relay stations in the service area of 
each main station. Continental transmitters have 
been allocated a single number from 201 onwards 



with each country having a block allocation giving 
room to add extra transmitters when required. United 
Kingdom transmitters are stored in film blocks 100- 
1444 and continental transmitters in blocks 1501 — 
3000. 

The data for each station are provided on the 
forms shown in Fig. 1 and, with the exception of 
the offsets, the details are assumed to be the same 
for all the channels at each station. The essential 
information is: 

(i) Transmitter number 

(ii) Transmitter location, i.e. latitude and longi- 
tude. The storage programme provides for the 
location to be specified as a National Grid 
reference, if required, and contains a routine 
to convert this to latitude and longitude by 
the method laid down in an Ordnance Survey 
booklet. 3 

(iii) Channel numbers ) Packed so that each 

> pair occupies one 
(iv) Corresponding offsets 'computer word. 

(v) Polarization 

(vi) Height of transmitting aerial above mean sea 
level 

(vii) Terrain clearance angles at transmitter at 
10° intervals of azimuth 'packed' so that five 
occupy one computer word. 

(viii) The horizontal radiation pattern of the trans- 
mitting aerial in terms of relative field strength 
at 20° intervals. This is converted to e.r.p. 
(dB relative to 1 kW) by the storage programme, 
and 'packed' so that six occupy one computer 
word. 

(ix) The numbers of up to fourteen test receiving 
sites (for United Kingdom transmitters only). 

All the above information is stored on one 
block of film (A). 

4.7. Test Location Data 

Each test location in the country has been 
allocated a number from to 4095. This defines 
the block number of film (B) on which the data is 
stored. 

The data consists of: 

(i) Location i.e. latitude and longitude but speci- 
fied as a National Grid reference if required 
(see Section 4.6.) 



(ii) Height of receiving aerial above mean sea level 
assuming always a receiving aerial 9-1 m (30ft) 
above ground level. 

(iii) Terrain clearance angles at 10° intervals of 
azimuth. 

4.8. Convergence 

The 'convergence* is the angle between the 
north/south grid line and true north at any point. 

4.9. Calculation of Distance and Bearings between 
two Locations 

The calculation of distance and bearings be- 
tween two locations is illustrated in Fig. 2(a). Let 
(#a^a) (#b^b) be the latitudes and longitudes of 
two points A and B, X be the bearing of B from A, 
Y be the bearing of A from B, and Z be the great 
circle distance AB. We have: 



tan Vi(Y + X) = -cot 1 /2(<6 B -<£ A ) x 



tan y 2 (y - X) = -cot V2(<f> B -<t> A ) x 



sintt(fl B -g A ) 

cosy 2 (0 B +<9 A ) 

cos%(# B -# A ) 
sin %(0 B +0 A ) 



sin y 2 (y-X) 

tany 2 Z = tanV 2 (tf B -0 A )x--^^ 
unless 6 B = # A wnen 



Z n cos 'A(y-X) 

tan— = cot V2(0 B + °a) x 

2 cos Vi(Y+X) 



If A = B then X = 0°, Y = 180° if 6 B > d k 
orX= 180°, y =0° if d B <d A 
andZ = \6 B ~9 A \ x 111-136 km 

The only places where these formulae can break 
down are: 

<% + ^A = °> '- e - A < B ' n different hemispheres 

or 6 Q + e k 

= 90°, i.e. A, B at North Pole 

2 

Neither of these cases is relevant to the U.K. 
problem. 

4.10. Calculation of Land and Sea Path Lengths 

The 'bearing' of the test location from the 
North Pole is first determined, this being defined 



as the angle west of 10°E. Its distance from the 
Pole is found simply since it is along a line of 
longitude, i.e. the distance is 111-136 x (90°- 
latitude) km.* The 'bearing' and distance of the 
unwanted transmitter from the North Pole are deter- 
mined in a similar manner. The bearings and 
distance between the unwanted transmitter and the 
test location are then calculated by the method 
described in Section 4.9. 

The method of calculating the land and sea 
paths is illustrated in Fig. 2(b). The great circle 
path between the transmitter and the test location 
crosses a family of coastal grid-lines of which, for 
clarity, only two are shown in the figure. The 
programme now goes through, in sequence, each grid 
line which falls between the unwanted transmitter 
and the test location, starting with that nearest the 
transmitter. If no grid line is crossed between 
these two points, then the nearest grid line is 
taken to represent the composition of the propaga- 
tion path. The distances along the line between 

(90 — latitude) is here expressed in degrees and 
decimals of a degree. 




transmitter 



Fig. 2 - Calculation of distance and bearings be- 
tween two locations 



each coastal intersection between transmitter and 
test location are multiplied in the ratio of the true 
great circle distance to that which is projected on 
the grid line. 

If the special case described above does not 
occur, the sequence illustrated by Fig. 2(b) is 
followed. The distance from the transmitter along 
path TR to the intersection with the first grid line 
(L) and the distance of L from the North Pole are 
calculated using the spherical triangle. If the tri- 
angle has angles A, B, C and the opposite sides 
are a, b, c expressed as angles subtended at the 
centre of the sphere, then: 

cot a. sin b = cos b. cos C + cot A. sin C 

In this case, side b, and the two angles A and C 
are known, so that a and similarly c are derived. 



calculation is now carried out for 1% time land, 
10% time land and, if the land-sea fraction is not 
zero, for 1% time sea and 10% time sea. 

The basic field strength Ff, is considered in 
terms of free-space loss minus an urban loss. We 
then allow 

(i) A loss Fi increasing linearly with the distance 
between transmitter and receiver, with an 
additional loss F 2 introduced at the horizon 
if they are screened from each other by the 
curvature of the earth. The F-l and F 2 losses 
depend on the path constitution, i.e. land or 
sea, and also on the percentage time for which 
the field is required (see Part 1 1 Section 3.2.) 

(ii) A loss (F ) depending on the terrain clearance 
angle within 15km of the terminals 



The process is repeated with the next grid 
line, resulting in distances c' and a'. The path 
between L and.L' is assumed to be represented by 
the portion of the grid line between distances c and 
c'. These land and sea dimensions are multiplied 
by (a-a)/(c-c) to give the correct distance LL". At 
the grid line immediately before the test location, 
the second point is the receiver instead of an inter- 
section point. The result of the process described 
above is a series of distances from the transmitter 
along the path TR defining the coastal intersections. 
It is arranged so that both transmitter and test 
location are on land, i.e. there is an even number 
of intersections. 



(iii) A gain (F ) depending on the terrain clearance 
angle and the aerial height 



The distance to the horizon can only be estima- 
ted approximately since mean land height is above 
sea level. For the purpose of this report, mean 
land height has been taken as 120 m above ordnance 
datum. This value of average land height is some- 
what arbitrary, but the precise value does not 
affect the calculation to an important extent. Where 
the intervening path is sea, an horizon height of 
zero is used. 



It was found when using this method that a 
false intersection at small intervals could be obtain- 
ed when the path is parallel to the coast-line. For 
this reason, it was decided to regard as sea astretch 
of land projecting into the sea but less than 5km 
wide, and vice versa. 



4.11. Field-Strength Calculation 

The information required for the field-strength 
calculation (i.e. distance between transmitter and 
receiver, height and terrain clearance angle at 
transmitter and receiver, and path constitution in 
terms of land or sea), is now available. 



The land-sea information is used to express 
the path constitution as a ratio of 'all sea', e.g. 
all land, 0; half and half, 0-5. The following 



The sum of the horizon ranges (D H ) is then 
given by 

D H = (2R) /2 i(h'j) y2 + (h' R ) A \ 



where h'j, A'r are the heights of transmitter and 
receiver above the datum height (always > 0) and R 
is the radio radius of the earth; R is given values 
depending on the land-sea path ratio and also on 
the percentage time for which the field is required. 



Because the terrain clearance angle takes into 
account details of the ground within 15km of the 
site, complete profile information is available for 
path lengths between transmitter and receiver of 
up to 30 km. 

Where the distance between transmitter and 
receiver is less than 30 km, or greater than 30 km 



8 



but within the horizon, the terrain clearance angles 
are redefined relative to the line joining trans- 
mitter to receiver, as shown in Figs. 3(a) and 3(6); 
these modified angles a x and a 2 are referred to as 
'terrain correction angles.' This means a correction 
to the terrain clearance angle of 0j degrees at the 
transmitter, where: 

T = id/2R + (h r - hR)/d\fl8Q/A 

where d is transmitter-receiver distance, and d R 
degrees at the receiver where: 

R =/d/2R +{h R -A T )/di|l80/ffJ 

The angles dj, 0r are those between the hori- 
zontal and the line-of-sight at the transmitter and 
receiver respectively. Thus, the terrain correction 
angle oti at the transmitter is given by 

Oft = {0 X -90) -0 T 

Similarly, at the receiver 

ct 2 = (0 2 -90) -0 R 



If the distance is greater than that to the horizon, 
as in Fig. 3(c), the terrain correction angles are 
those relative to the line from aerial to the horizon. 
The correction angle 0j degrees is given in this 
case by: 



T = (2h r /R) 



|l80/77> 



In order to make some allowance for ground 
clearance, the terrain correction angle has to be 
greater than 0-1326° for the path to be considered 
clear. This value was selected after investigating 
measurements and related factors, but is also 
consistent with first Fresnel zone clearance at 
the mid-point of a 30 km path for a frequency of 
about 1000 MHz. 

If the path between transmitter and receiver is 
less than 15km, and either a ± or a 2 is greater than 
0-1326° then there must be line-of-sight. The field 
strength is then F b as defined above. 

If neither terrain correction angle represents a 
clearance, then the linear loss F a is applied, 
together with the loss F shown in Fig. 4 defined 
by an angle a defined as fol lows: 

a = ifa 1 >0 and a 2 >0 

a = ati ifai<0 and o< 2 >0 

a = cc 2 ifai>0 and a 2 <0 

a = a>. + a 2 ifa 1 <0 and cc 2 < 

If the path length is between 15 km and 30 km, 
and both terrain correction angles represent a 



clearance then the field strength is the basic 
figure Fb- " either terrain correction angie is less 
than 0-1326° then the linear loss F x and the obstacle 
loss F for the angle equal to the sum of the 
negative terrain correction angles are applied to the 
basic field strength F b . 

If the path length is greater than 30 km, then 
the linear loss Fi is applied together with the 
additional loss F 2 if transmitter and receiver are 
beyond the horizon from each other. At each 
terminal, if the terrain correction angle is less 
than 0-1326° then the terrain correction loss (F ) 
is read off from Fig. 4. 

Terrain correction gain is applied only to path 
lengths exceeding 30km and is determined by the 
angles shown in Fig. 3(d). From each terminal we 
draw a line to the horizon if the path length equals 
or exceeds Dh> or to the other terminal if the path 
length is less than Dh. We also draw a line from 
the terminal to a point at the horizon height, at 
15km from the terminal; och is then defined as 
either the angle between these lines, or as 0-1326° 
if oc H is less than 0-1326°. 



horizontal — 




horizon height (taken as 
120 rn tor land, Om for sea) 




horizon height (taken 
as 120m for land, Om for sea) 
Id) 



15 30 SO 

path length, km 



Fig. 3 - Path geometry diagrams 

(a) Terrain correction loss; path length <30km (aj, a 2 

are negative in this figure) 
(fa) Terrain correction loss; path length >30km <D// 

(c) Terrain correction loss; path length > D// 

(d) Terrain correction gain; derivation of angles and 
interpolation (one terminal). (<Zi, a.f] are positive in 
this figure) 

a \, 02 are terrain clearance angles 
a l, a 2 are terrain correction angles 
Df/ is the sum of horizon distances 



For path lengths between 30km and 60km the 
angle used is a value obtained by a linear inter- 
polation, as indicated in Fig. 3(d), between the 
terrain correction angle oi (or 0-1326° if a t is less 
than 0-1326°) and the value of on determined as 
described above. 
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assuming a coefficient of correlation (p) of 0-5. 
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Fig. 4 - Terrain correction 

The terrain correction gain is obtained from 
Fig. 4 for each end of the path, and the gains added. 

Having calculated 1% time and 10% time field 
strengths for land, and if necessary for sea, the 1% 
time and 10% time figures for the appropriate mixed 
path are derived as proportions between all-land 
and all-sea. 

4.12. Calculation of Median Field Strength ando- T 

Let the two median- location field strengths 
exceeded for 1% and 10% time be Ei >50 and E i0 50 
respectively. Assuming a Normal distribution in 
time, the median value (M), and the standard devi- 
ation (ctt) in time are calculated using equations:- 

cr T = 0-9559 (E 1>5 o-E 10 ,5o) 
M = 2-2259 Eio lS o- 1-2259 E lj50 

The a j for the interfering signal level relative 
to the wanted signal level is calculated from the 
individual values ctju and orw from the formula: - 

a T = S^TW - 2 P CT TU<7TW + a lU> 
assuming a coefficient of correlation (p) of 0-5. 

Similarly a'\_ for the interfering signal level 
relative to the wanted signal level is calculated 
from 



It is also assumed that o\_ for any signal is 5dB 
hence a' L = 5dB (See Part 2 2 , section 4.3.). 

4.13. Calculation of E.R.P. Factor 

The e.r.p. term for the unwanted station in the 
direction of the test location is found using linear 
interpolation between the 20° values stored after 
'unpacking.' 

4.14. Carrier Frequency Offset Protection Ratio 

Apart from nominal zero offsets the only offset 
conditions considered are ±2/3 line-scan frequency. 

The protection ratio required when there is an 
offset between wanted and unwanted signals is 
30dB. 

The protection ratio required for no offset 
between wanted and unwanted signals is 45 dB. 

4.15. Calculation of Receiving-Aerial Directivity 
Factor 

The receiving-aerial directivity factor is a 
function of the difference in bearings of wanted and 
unwanted transmitters from the test location. The 
polarization of each transmitter is already available 
in the store, and an allowance is included in accord- 
ance with the information contained in Fig. 5. 
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Fig. 5 - UHF receiving aerial discrimination curves 
Angle relative to direction of main response 
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4.16. Median Protected Field Strength, M' 
The median protected field strength (M) 
(median field strength, M) + (e.r.p. term) 
+ (protection ratio for carrier frequency offset) 
+ (receiving aerial directivity factor) 



4.17. Extra Protection Ratio for Long-Term Inter- 
ference 

The 95% time protected field strength is 
calculated from the formula 

E '95,50 = M '+ 1 -645<7x 

If (E 9 s,50 -MO < 10 (i.e. ct t < 6.1) then the 
particular protected field strength is increased by 
the extra protection ratio (10 - 1-645a T ) dB. 

4.18. Multiple Interference Calculation 

Since the two variables Involved, time and 
location, are independent, the standard deviation 
a TL of tne f^'d strength for a combined variable 
time-location, is given by 

° TL = a T +° L 



The combined effect of many interfering sources 
is found in the following manner, with the object 
of determining the median field strength from the 
wanted station at the test location that gives a 
protection probability product in time-location of 
95% ± 1%. An initial field strength of 70dB is 
assumed and the probability of protection against 
each interfering source acting separately is found. 
The product of these probabilities is calculated 
giving the probability of protection with this field 
strength against all sources. In general this will 
not be the value of 95% ultimately required. A note 
is made of the a' T! _ of that source which affects 
the product to the greatest extent, i.e. the source 
whose probability of protection is least. The 
assumption is made that this a 'j\_ is representative 
of the slope of the probability product curve at this 
field strength. It is then used to determine from 
the probability product for 70 dB the field strength 
for 95% probability of protection. This new estimate 
of the field strength required is now used to find a 
new probability product and the process repeated 
until a field strength is found which gives a prob- 
ability product in time-location within 1% of the 



required value. If, when the third or subsequent 
value of field strength is found, the field-strength 
sequence is oscillating, i.e. E\n) - E'(n-1) has a 
different sign from E'(n-1) - E\n-2), a speeding up 
of the convergence is obtained if the next value of 
field strength used is VilEln) + Efa-I)] instead of 
E\n). If after six loops, the test is still not satis- 
fied, then a linear interpolation based on the dif- 
ferences between the achieved and wanted prob- 
ability products for two loops is used to derive the 
next value of field strength. 

Having determined E' 9S , the whole process is 
repeated for 90%, 82%, 70% and 56% probability 
products, giving E 95 , E 90s E' 82 , E' 70 , E 5S . Each 
interval e.g. E95 to E 90 , is considered in turn. 
The bounding field strengths are assumed to define 
a Normal distribution and since the probability 
products at these field strengths are known, the 
median M' and standard deviation a'j{_ for this 
distribution are found. 



The mean field strength in the interval is taken 
and the probability of protection against each inter- 
fering source, when the variation of field strength 
with location only is considered, is found. The 
product of these probabilities is calculated and, 
together with the previously calculated median field 
strength, is used to derive a'\_ for the interval. This 
i?'l is limited to a maximum of 5dB. Hence, from 
o-'-fL and ct' l , o-'-r is calculated. From the values 
of median and a'y, the field strength required to 
protect against interference for 95% time and 50% 
location is found. This is checked to see if it lies 
between those values of field strength used to 
define the interval. If it does not, the next interval 
is considered and so on. 

Having found the interval such that the field 
strength required to protect for 95% time and 50% 
location lies within it, the field strength required 
to protect for 95% time and 70% location is calcula- 
ted using the appropriate a'\_ and assuming this may 
still be referred to the parameters of the same 
interval. 

The other information required is the percent- 
ages of time for which field strengths of 70 dB and 
80 dB protect 70% of locations against interference. 
Again the interval into which the field strength 
falls is found, and the appropriate a'j are used to 
calculate the percentage of time. 

A typical computer print out is shown in Table 
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TABLE 2 
Co-Channel Interference Calculation with Terrain Correction 



annel 


Date 




Service Area 


26 


2/9/66 




C P 


Test Location 


95% Time 


95% Time 


70% Locn 


70% Locn 


Dist. Bearing 


50% Locn 


70% Locn 


70 dB 


80 dB 


(km) (Deg) 


F.S. (dB) 


F.S. (dB) 


% Time 


% Time 


From TX 


Prot. 


Prot. 


Prot. 


Prot. 



224 



35-6 



191-5 



61-0 



63-3 



98-6 



99-9 



TX 



Stratfd 

Ramsbot 

Skipton 

Halifax 

Swansea 

Rhondda 

C Risca 

Hexham 

Cromer 

Alburg 

N Yorks 

Wrekin 

Stockld 

Ayrshir 

Banff 

Lewis 

Ireland 
Kippure 
Mulrany 
S Naght 

Germany 
Bonn 
Fulda 
Minden 
Sleswig 
Stutgat 

France 
Besacon 
Bourges 
Meziere 
Nancy 
Nantes 
Rouen 

Holland 
Lelstad 

Belgium 
Molhout 

Norway 
Flo 
Oevre 



50% Time 

50% Locn 

F.S. (dB) 

Prot. 

14 

-14 

-33 

-12 

-23 

-10 

21 

-55 

13 

34 

28 

20 

36 

-60 

-91 

-134 



RELATIVE 
S.D. 
Time 



-12 
-88 
-86 

-61 

-105 

-74 

-85 

-113 

-87 
-35 

-3 
-58 
-18 

25 

1 

-28 

-184 
-132 



6 
13 
15 
13 
12 
10 

8 
20 
10 

7 
15 

9 

8 
26 
36 
42 

25 
36 
33 

25 
34 
33 
50 
35 

31 
24 
18 
26 
27 
14 

26 

20 

73 
63 



CONTRIBUTIONS 
95% Time 
50% Locn 
F.S. (dB) 
Prot. 
24-3 
8-3 
-7-8 
10-5 
-2-7 
6-9 
35-9 
-21-7 
29-8 
47-3 
54-0 
36-1 
49-9 
-17-2 
-31-9 
-64-6 



30 
-28 
-31' 

-19. 
-49 • 
-19. 
-2 
-54> 



-36-2 

5-1 

27-7 

-15-4 
25-9 
48-9 

44-9 

5-2 

-62-9 
-26-9 



Aerial 
direct. 
(Deg) 

303 
321 
325 
327 
271 
273 
273 
333 
13 
32 
337 
305 
249 
32 
334 
322 

288 
287 
304 

79 
79 
62 
44 
95 

118 
140 
99 
106 
179 
143 

53 

73 

4 
16 



TX 

Bearing 

(Deg) 

314-6 
332-9 
336-8 
338-0 
282-6 
284-2 
284-8 
344-2 
24-2 
43-5 
348-8 
316-4 
260-3 
331-5 
345-7 
333-5 

299-5 
298-6 
315-0 

90-9 
90-0 
73-6 
55-2 
106-1 

129-8 
151-9 
110-5 
117-3 
190-8 
154-5 

64-8 

84-1 

15-8 
27-0 
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